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O ■ Abstract 



We consider a system of N weakly interacting fermions with a real analytic pair interaction. 
We prove that for a general class of initial data there exists a fixed time T such that the difference 



between the one particle density matrix of this system and the solution of the non-linear Hartree 



equation is of order N 1 for any time t < T. 
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1 Introduction 

The Hartree-Fock theory is a fundamental tool in atomic physics, chemistry, plasma physics and 
many areas of quantum physics. It is also an important numerical instrument to calculate atomic 
and molecular structures. Despite numerous applications of the Hartree-Fock theory, many basic 
theoretical questions remain unsolved. One area where significant progress was made concerns the 
ground state energy of large atoms and molecules. Consider the simple case of a neutral atom with 



*Partially supported by NSF grant DMS-0200235. On leave from School of Mathematics, GeorgiaTech, USA 
^Partially supported by NSF grant DMS-0307295 and MacArthur Fellowship. On leave from Courant Institute, 
New York University, USA 



1 



nuclear charge Z. It was first proved by Lieb and Simon jlUl lllj that the Hartree-Fock theory gives 
the correct asymptotic energy to the leading order Z 7 / 3 as Z — > oo. The next important step came 
more than a decade later as Bach [I] proved that the error between the Hartree-Fock and the true 
atomic energy is less than Z 5 / 3- " 5 for some small 5 > 0. Similar result with a very different method 
was announced in jlj and was proved in jHJ. 

The goal of this paper is to justify a time-dependent mean-field theory for the evolution of 
interacting fermions under a weak pair interaction with initial data localized in a cube of size of 
order one. The last restriction actually provides the length scale of the system. The interaction 
potential varies on the same length scale. While one might want to add a background potential, 
we shall keep the model simple to focus on the many-body interaction effect. We work in d = 3 
dimensions, but our result holds in any dimension. The Hamiltonian describing such a system is 
given by 

2 N i 

3=1 j,k 

acting on /\^ V L 2 (M 3 ), and the Schrodinger equation is given by 

iedti> t = H N ii t . (1.2) 

Here we have chosen the strength of the interaction between fermions to be of order 1/N. Examples 
of such systems with a small coupling constant can be found in astrophysics and plasma physics. 
For gravitating systems, the strength of the interaction is dictated by the gravitational constant 
and thus the mean field approximation is suitable. The Coulomb singularity, however, is difficult to 
control. If one wishes to use (jl.lj) to model the dynamics of white dwarfs, the kinetic energy has to 
be further modified to be the relativistic one, according to the famous observation by Chandrasekhar 
[H], see a rigorous account in [Ej. For the plasma physics application, the weak pair potential models 
combined electron-electron and electron-background interactions. 

From now on we will fix a particular relation between e and N, which is motivated by the 
following argument. We model a system of N fermions at energy comparable with the ground state 
energy of the system. The potential energy per particle is of order one. It is well-known that the 
kinetic energy per particle of N fermions, i.e., — ^e 2 A Xj , in a cube of size one scales like e 2 N 2 / 3 in 
the ground state. In order to keep the kinetic energy per particle of order one, we need to choose 
e = iV -1 / 3 , a convention we shall use for the rest of this paper. With this choice, the kinetic and the 
potential energy per particle in Hn are comparable. This is the basic physical criterion to obtain a 
limiting dynamics (as N — > oo) that captures the nonlinear effect of the interaction. Notice that we 
have kept the free evolution in the form of iedtip = — ^e 2 A^ so that the free evolution has a limit as 
e — > 0. The equation (jl.2j) is formally semiclassical with a mean-field interaction potential at high 
density. Our choice of scaling is the same as in ^1] and (The interpretation of the origin of this 
scaling in Jl] is somewhat different.) 
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In order to take the limit e — > 0, we need to recast the Schrodinger equation using the density 
matrix. For any wave function ipN,t, define the corresponding density matrix by 7jv,t = 7t^p Nt , where 
7r^ = \ip)(vp\ is the orthogonal projection onto ip. The kernel of 7^ is then given by 



7A r ,t( x >y) = i>N,t{*)il)N,t{y) ■ (i-3) 

The notation x typically stands for x = (x±, ... ,xn). Depending on the context sometimes it may 
denote a shorter vector of x's. 

We recall that a self adjoint operator 7 is called density matrix if < 7 < 1 and Tr 7 = 1. If 
the density matrix of the system is a one-dimensional projection then we say the system is in a pure 
state, otherwise it is in a mixed state. The Schrodinger equation (|1.2|) is equivalent to the Heisenberg 
equation for the density matrix: 

ied t j N)t = [H N , lN>t ] , [A, B] = AB — BA . (1.4) 

(n) 

The n-particle density matrix, j Nt , is defined through its kernel 

7nI(xi, -.,x n ;yi, ..,y n ) := J dx n+1 ..dx N jN,t(xi, .., x n , x n+ i, .., x N ; yi, .., y n , x n+1 , .., x N ) (1.5) 

(n) 

for 1 < n < N, and r y N t := otherwise. Define the Wigner transform of the one particle density 
matrix in the scale e by 

W ( m\x;v):=j^ J e-^^x + e^x-e^jdrj. (1.6) 

Recall the nonlinear Vlasov equation for a phase space density /: 

dtft{x,v)+vV x f t {x,v) =V x (U*Q t )-V v ft(x,v) , (1.7) 

where 

Qtix) := J f t (x,v)dv 

is the configuration space density. It was proved by Narnhofer and Sewell |14j that WpP converges 
weakly to a solution of the Vlasov equation (|1.7|) provided that the Fourier transform of the potential 
is compactly supported, in particular U is real analytic. The regularity assumption was substantially 
relaxed by Spohn [T5] . 

Define the Hartree equation for the time dependent one-particle density matrix uit by 



ied t Ldt 



e 2 



—A + U* 6t ,uj t (1.8) 
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where Qt(x) '■= Ut(x,x). Note that the Vlasov equation Q1.7JI is the semiclassical approximation of 
(|1.8|) . One can extend this equation to the Hartree-Fock equation by including the exchange term 

-2 



iedtUJt 



U(x- z) — U(y- z) u t (x, z)u t (z, y)dz . (l.l 



Our main result proves that the Hartree equation correctly describes the evolution of the Schrodin- 
ger equation (|1.4jl up to order 0(e). More precisely, it states that for short semiclassical time the 
difference between the Wigner transform wjp t of the solution to the Schrodinger equation (|1.4|) and 
the Wigner transform of the solution of the Hartree equation (jl.8|) is of order 0(e 3 ) provided that 
the potential U is real analytic. In other words, all e 2 corrections come from the difference between 
the Vlasov equation (|1.7|) and the Hartree equation (|l.Hjl : hence they are related to the accuracy of 
the semiclassical approximation in the one-body theory. In particular we show that all correlation 
effects are of order at most 0(e 3 ). 

In fact, the main correlation effect, the exchange term, is expected to be order e 3 for smooth po- 
tential and e 2 for the Coulomb potential. Our interpretation of the Hartree-Fock equation resembles 
the theory concerning the ground state energy for atoms where the Hartree-Fock theory is proved 
to be correct up to e 2+s smaller than the leading term pQ. The analyticity condition and the short 
time restriction of our result is unsatisfactory; it nevertheless shows what the correct formulation of 
the time-dependent Hartree and Hartree-Fock theories should be. 

In order to see the effects of the exchange term, i.e. to show that (|1.9|) approximates the quantum 
dynamics even better than (|1.8j) . we would need to consider e 3 correction for the smooth case or the 
e 2 correction for the Coulomb potential. Notice that our approach is perturbative and in principle all 
e 3 corrections, including the exchange terms, can be calculated. However, there are other sources of 
e 3 corrections (see the last two terms in Q4.15JI in Section 4) which make the exchange correction less 
prominent. This should be compared with the Coulomb case where all e 2 corrections are expected 
to be from the semiclassical approximation to the Hartree equation and the exchange terms. 

In a recent paper Gram et al. [Hj proved the convergence of the Wigner transform wffl (of the 
solution of the Heisenberg equation Q1-4JI ) to the solution of the Vlasov equation under the assumption 
that the initial wave function is of the semiclassical form tp = Ae %s ^ e . The result also provided 
error control and the proof is carried out by concise inequalities as opposed to weak convergence 
method in ^1] and |15j . The main restriction is the initial wave functions to be of the semiclassical 
form. Although this type of wave functions is suitable for bosons, fermionic wave functions are 
antisymmetric and thus vanish frequently. Notice that in the neighborhood of the zero set of the 
wave functions, the semiclassical approximation is difficult to apply. In particular, one naive attempt 
(for fermionic case) is to choose S symmetric and A antisymmetric. But J |VA| 2 will be of order 
iV 3 / 5 and this violates a key assumption in this paper. 

The recent work of Bardos et al. 2 considers the equation 



(-a^A x . + ^J2 U ( x j- x k))^N,t (1-10) 

V 3=1 j,k / 
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with an arbitrary a > (we have put h = 1 which is a constant of order one in this paper). In the limit 
N — > oo, it was proved that the difference between the one-particle density matrix 7^ = |^>;v,i)(V'iV,t| 
and the solution to the corresponding time-dependent Hartree-Fock equation vanishes in the trace 
norm provided that the initial data is a Slater determinant (and some other assumptions). Notice 
that the time scale in (|1.1U|) is of order e = iV -1 / 3 smaller than (|1.2|) . Thus for initial data considered 
in ^1] ^S] and the present article, the one particle dynamics of (jl.lU|) is governed by a free evolution; 
the effect of the interaction given by U vanishes in the limit N — > 00. We shall make a more detailed 
comparison in Section 3. 

Finally we comment on the method. Our approach is to based on the BBGKY hierarchy and 
iteration scheme. There are two major elements in the proof. The first one is the control of error 
term. Since we work on the BBGKY hierarchy for finite N, we need to control the error term in the 
iteration scheme. Here we used that the trace norm of the density matrix is preserved. The second 
observation concerns the combinatorics. As usual, the BBGKY hierarchy will produce a n! factor 
under iteration. However, in the setting of this paper, there are extra sources of n!, for example, we 
will need to take high moments of the interaction: 

J \U(Om m d^C m m\ (1.11) 

See 1)1.11(1 for precise assumption. Since time ordered integration provides only a l/n\, we will have 
to prove that the combined effects of the factorials from all sources is just a single n\. See the proof 
of Lemma 4.1 for details. 

Acknowledgments. We are grateful to Herbert Spohn for useful discussions. 



2 Notations 

We first fix the notations and recall some definitions. The n-particle density matrix 7^ is defined 
through the equation (|1.5|) and clearly satisfies the following normalization 

It is well-known that the one particle density matrix satisfies the following operator inequality (see 

M) 

< 7#J < ^ • (2-2) 



Therefore, we can write 7^ as 



1 00 
(i) 1 \^ 



N 
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where Hj is the orthogonal projection onto tpj, and where a,j £ [0, 1] for all j £ N with Y^jLi a j = 1- 
Note that in the definition of the re-particle density matrices we followed the convention that the trace 
of the density matrices is normalized. In standard iV-body theory an additional N(N — 1) . . . (N — 
n + 1) factor would be present in (jl.5[) . 

2.1 Wigner Transform 

The Wigner transform of an iV-body density matrix 77\r(x;y) is defined by 

Wjv(x;v) := q^ 3N J e _JV ' y 7Jv(x + |,x- |)dy . (2.3) 
From Tr 77V = 1 it follows that 

J dxdvw;7v(x,v) = 1. (2.4) 

Since the velocities of the N particles are of order iV 1 / 3 = e~ 1 , we rescale the Wigner transform wn 
so that its arguments be typically of order one. Thus we defined the rescaled Wigner transform by 

^, e (x,v)=Ty w (x,v):=e- 3 ^(x,v/e) = ^-L F J dy 7 iv(x + - ^) e -^. (2.5) 

The factor e~ 3N guarantees that the normalization Q2.4JI holds for the rescaled Wigner transform 
Wn,e( x , v ) as well. The inverse transform is given by 

77v(x,y)= J W N , £ (^^,u)e^-y> u / £ du. 

In particular, the particle density at the point x is given by 

p(x) := 7at(x,x) = J Wjv )£ (x,u)du. 

In this paper, we are concerned with the rescaled Wigner transform only, so we shall drop the 
adjective "rescaled" and the e index from the notation. The rescaling parameter e will always be 
related to the total number of particles as e = iV -1 / 3 . 

The time evolution of the Wigner transform Wjv(x, v) is given by the Wigner equation 

N 

d t W N {t; x, v) + Vj ■ V X] W N (t; x, v) 
i=i 

e iy - (u - v) W iV (t;x,u)dudy, (2.6) 



IE' 



(2VT 



,37V 



U x + 



U 



£y 
2 
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which can be easily derived from the Heisenberg equation (|1.4|> . 

It is tempting to consider the Wigner transform as a probability density on the phase space. 
The problem with this interpretation is that Wn(x, v) is not positive. In order to make the Wigner 
transform positive we may take convolutions with Gaussian distributions. We can define the Husimi 
function by 

where * x denotes the convolution in x-space and 

is the centered Gaussian distribution in 3iV dimensions with variance 5. It is easy to check that 
-Hjv 52 — ^ ^1^2 > £• The Husimi function is normalized according to 



J F^' 52 (x,v)dxdv 



Ml 



and thus can be considered as a probability density on the phase space. The accuracy of the H% ' h > 
is of order 5\ for the space variables and 82 for the velocity variables in semiclassical units. 

As a side remark, we recall that for 5\ := S, 62 ■= £<5 _1 the Husimi function is just the standard 
Gaussian coherent state at scale 5: 

H^ 6 ~\x, v) = CW,i(x,v) := (27re)- 3A ^7V,<>iv) 

where tt^ v = |0xv)(^xvl * s t ne orthogonal projection onto the state 



(7T(5 2 ) 3iV / 4 25 2 

(k) 

The fc-particle Wigner transform Wj^ is defined to be the Wigner transform of the k particle 

(k) 

density matrix 7^ . Clearly, it can be viewed as the A;-particle marginal of Wn since it satisfies 
W^\xi,...x k ;vi,...v k ) = J J dy 7 JJ } (x+y;x- e~ JV ' y 

(2.7) 

Wjv(xi, . . . x N ;vi, . . . v N )dx k +i ■ ■ ■ dx N dv k +i . . . d-u^ 



The are normalized as 



y W r Jf ) (x,v)dxdv = 1. 



Notice this definition is consistent with the definition (J1.6p • We now give some examples of ./V particle 
wave functions. 
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2.2 Some Examples 

(n) 

One of the most important assumptions of our results is that the n-particle Wigner transform W N of 
the initial state is factorized in the limit N — ► oo. At first glance this assertion might seem surprising, 
since we are dealing with a system of fermions. In the following we present some typical situation 
where this condition is indeed fulfilled and we show a very atypical example where factorization is 
wrong. 

The standard examples of many-body fermionic states are the Slater determinants and the 
quasifree states. 

1. Slater determinants. For any orthonormal family {<~fj,j = 1, . . . , N} define the determinant 
wave function 

N „ N 



j=i ViV -aes JV j=l 

The one particle density matrix is given by 

1 " 

7 (1) (x,x / ) = —^2ip j (x)ip j (x') . 

3=1 

The two-particle density matrix is 

TV 



7 (2) (x,y,x',y / ) = 2jy ^ _ ^ ^ [<Pk(x)<pe(v) ~ t Pk(v)<Pt( x ) VkWMl/) ~ <Pk(y , )<Pt(. x ' 

(x, x')7 (1) (y, y) - 7 (1) 0, y')7 (1) (y, z') 



(2.8) 



iV- 1 

£. Quasifree states. An A^-particle state w is called quasifree, if its fc-particle density matrices 
factorize by Wick theorem 

J k \x 1 ,...,x k ;y 1 ,...,y k ) = j— _ - — _ - - ^ det(wW(Xj,yj))j=i,...,fc • 

The unusual prefactor is present due to our choice of normalization. In particular, quasifree states are 
characterized by their one particle marginals. For example, Slater determinants are pure quasifree 
states. 

The concept of quasifree state can be generalized to grand canonical states with an indefinite 
particle number. In particular, any normalized density matrix 7 on L (R 3 ), can be realized as 
a fermionic quantum state whose one particle density matrix is 7. The state can have expected 
particle number up to 1/||7||. In fact, with N := J dxTr (uja x a x ) we have j(x,y) = N~ lr Fr {ua x a y ), 
where at, a x are fermionic creation and annihilation operators. Thus, for any tjj G L 2 (IR 3 ), 

WM) = ^Tr (ua^) = -llV (ua^a\) + M! < Ml! (2 . 9) 
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where the operators and annihilate and, respectively, create a fermion in the state ip. 

Example 1: Let := [0,2-7r] 3 and consider the states fk{x) = (2ir)~ 3 / 2 e lkx x{x € f2) with 
\k\ < cN 1 / 3 , k G Z 3 . The number of states is 0(N). We consider the pure state \P = /\ of the 
./V particle system and we compute the marginals of its Wigner transform. The one particle density 
matrix is given by 

7 (1) (x,x') = ^ <Pk{x)<fk{x') 

k-.^KcN 1 / 3 ^ 

= ^ X(g ; 2 t ) l" ) E ')~ x( x,x' G ^)/((x-x')^/3) 

^ 11 ' k:\k\<cN^I^ 

with some decaying function /, such that /(0) = 1. Its Wigner transform with rescaling parameter 
e = JV-V3 i s 

= lUf E / dy e ^ e "^ ^ 72^3 x(x e n)x(M < C ) 

v 7 fcilfcl^cAn/s- 7 v ; 

when TV — > 00. Using (|2.8|) the two particle density matrix can be computed as well. Its Wigner 
transform is given by 

^ (SF^y / * — K> (* + f - f ) 7« (* + f - f 

«(, 1 + f„ 2 -f) 7 «( I2 + f„ 1 -f)} 

(2.11) 



7" 



Notice that the first term is W^\x\, v\)W^ l \x2, V2) after neglecting the error N/(N — 1) = 1. 

The second term is the so called exchange term and it vanishes as N — > 00. By H2.1U|) this term 
can be written as 

Wj^(x,v) = X ( Xl,X ^ e ^ / dy e -iv-y e ik(x 1 -X2)+ike(yi+y 2 ) e ie(x2-x 1 )+iis(yi+y 2 ) (2.12) 

(27r) iV fc ^ J 

Thus, for an arbitrary function J(x, v) = Ji(x)J2(v) we have 
y dxdvJ(x,v)W e ( ^(x,v) 

jV 

= ^pE ydxJ 1 (x 1 ,x 2 )e l ( fc - £ )^-^ydyJ 2 ( 2/1 ,y 2 ) e W)^ + ^) (2.13) 
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For any smooth functions J%, the x integration is very small unless £ ~ k. Thus the order of the 
exchange term is 1/N. Notice that if we take J(x,y) ~ \x — then the exchange term becomes 
of order iV~ 2 / 3 , consistent with standard pictures from semiclassical limits of atomic and molecular 
energies. Indeed, since 



dx e lkx -^— ~ — !— 
\x\ |fc| 2 



we obtain that in this case 



cN 

W ( 2 )(xv)~— V - N- 2 / 3 

lfe-€l=l v ; 

Instead of choosing ^ l \x,x') as in (|2.1U|) . we could also define 

7 (1) (x,x') = ^Y1 f(k)<f k (x)^) (2.14) 

for an arbitrary distribution f(k) with < f(k) < 1 for all k £,Z 3 , and with ^2^f{k) = N so that 
7^^(x,x') is a density matrix satisfying the conditions (|2.1j) and (|2.2|) . Typical distributions are of 
the form f(k) = g{ek) for some smooth function g. In this case the one particle density matrix is 
supported within a distance of order e from the diagonal, and the exchange term, analogously to 
()2.12j) . vanishes in the weak limit N —* oo. 

Example 2: Let w be a smooth decaying function. We define 

<p k (x) = e A '\ 



e 

where k runs over lattice sites with \k\ < c/e, and e = iV _1//3 , as always. In other words, (fk(x) 
represents a state localized inside a sphere of radius e around the lattice site k. It is a straight-forward 
exercise to show that the exchange term in the two-particle Wigner transform of /\ iff. is again of 
order 1/N . Indeed, we obtain in this case that 



f 1 N f —k 

/ dxJ(x)W c ( x 2 )(x, v) = —-J2 dxdyJ(x 1 ,^ 2 )e- iv -^(^— + 



xi - k yt\ fx 2 -k y 2 

UJ 



N 



2 ) V e 2 



'x 2 - £ y 2 \ (x\ -i yi 
x uj UJ 



2J \ e 2 



(2^) 6 iV 2 
v ; k,e=i 



j>~] J dxdyJ(exi,ex 2 )e lv ' y u(xi ~ ~ + y) ^ ( x 2 - ^ - y 



, i , V2\ ( £ yi 



(2.15) 
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Clearly, for any smooth function J(x), only terms with k = i give a considerable contribution to the 
sum; therefore the right hand side of the above equation is bounded by 1/N . 

Example 3: In the last two examples the kernel of the one particle density matrix, j^'(x,y), is 
concentrated on the diagonal \x — y\ < iV -1 / 3 . Suppose now that we are given a one particle density 
matrix ~f(x,y) on [0,27r] 3 x [0, 2ir] s , which satisfies < 7 < (const. )/N, Tr 7 = 1 and it is supported 
near the diagonal. Let 

%x, y) := (3 j(x, y) + -y(x + e, y) + j(x, y + e) + j(x + e, y + e) (2.16) 

with e := (0, 0, 2ir). We can choose the constant (3 so that Tr 7 = 1 and we still have 

< 7 < (const. )/N . 

Thus we constructed a density matrix which is not concentrated on the diagonal. The corresponding 
quasifree state can be constructed by standard procedure. 

This construction can be carried out on the level of the wave functions as well. Let tp k , k £ Z 3 , 
|^| < cN 1 / 3 be N orthonormal one body wavefunctions supported the cube [0, 2-7r] 3 as in Example 1. 
Then 7 := jj Yl,k is supported near the diagonal. Define 

i> = /\ipk , with ip k (x) := 2~ 1/2 {ip k (x) + ip k (x + e)} . 
k 

Then the one particle density matrix is of the form 7 from (|2.16|) . It is concetrated around three 
submanifolds x = y and x = y ± e and not just along the diagonal |x — y\ < iV -1 / 3 . In particular, 
the exchange term is still of order 1/N. 



The fact that the exchange terms are of order 1/N in all these three examples tells us that there 

(2) 

is a large class of initial data for which is factorized in the weak limit N — > 00. Similar result 
can be obtained for any n-particle function if n is fixed: 



lira Wp(xi,...,x n ,vi,...,v n ) = TT Wi'(xj,Vj) , (2.17) 



or, more precisely, | (J, Wjf 1 — {W^)® n )\ < c/N, where the constant depends on n and on the smooth 
test function J(x, v). 

2.3 The Fourier Transform of Wn(x., v) 

Instead of working directly with the Wigner function Wjv it is often more convenient to work with 
its Fourier transform, which we define as 
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where x and p = — iV x are the position and momentum operators on L 2 (M. 3N ). Narnhofer and Sewell 
defined the same quantity with a somewhat different notation (see (3.2) in |14j): the £,,rj variables 
are interchanged and the conjugate is considered. Noting that 



and hence 



(e-^-P+t-*^) (x) = e^e-^ipi* - ery), 

we have 

Mtv) = J e^ x 7 Jv(x-^,x + ^)dx = J dxdvW^(x,v) e -^ x -^ v 

^v(x,v) = J M^K^+^d^dt 

Notice the operator norm of e _ ^ £T, 'P + ^' x ^ is equal one. Since 77V is positive and Tr 77V = 1, we have 

\HN&ri)\ = \Tr lN e<^^\ < 1. (2.18) 

For any J(x, v) with ||./||L 1 (d£,dT7) bounded, we have 

\(J,W N )\ = \(J, m )\ < ||J|| Ll (d€,dT,), (2-19) 



where (/, g) := J dxdv/(x, v)g(pc, v). Therefore, one can always extract weak limit points of Wigner 
transforms. 

The time evolution of //jv(£, ij) can be easily derived from ()2.6|) : 

JV 

3=1 

N 

~ 2e2S ^Z I dq ^io) sin (f ~ ^(*' £i 3 •••> £j -?>•••> & + g, 6v; v), 

j<k J 

(2.20) 

where we defined U(q) := (27r) -3 f dxe~ iqx U(x). We denote by fj,£ the Fourier transform of the k 

(k) 

particle Wigner function . Then we have 

4>(ei ...6.://:. ..'/„) / dxdv < fc) (x,v) e - ix ^- iv -" 2 

= MJv(Ci 3 • • • > 6c, 0, . . . , 0; 771, . . . rj k , 0, . . . , 0), 
if k < N and [J$ = otherwise. From (|2.18l) 

\$\t,n)\<l (2-22) 

is valid for all k. 
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3 The BBGKY Hierarchy and the Main Result 

The family of marg inals {W^} n=i,..,N satisfies a hierarchy of equations, usually called the BBGKY 
hierarchy, which can be derived from (|2.20|) (also using the symmetry of wffl): 



X\ , . . . X n , V\ , . . . V n ) 



(n). 



3=1 



(2vr) 3 ™ i 



E /^^ + f-^-f)-^-f 



<j<k<n 



x e ! E^i wK' xi, . . . , x n , ui, . . . ,u n )dui&yi . . . du n dy, 

i(l — ne 



(2k) 



3n 



3=1 



^1 x i ^ ~ x «+i 



e *E "=!%•(%— «j0 



)duidyi . . . du n dy n du n+1 dx n+1 . 



(3.1) 



The main goal of this paper is to compare solutions of this hierarchy of equation with tensor products 
of solutions of the one particle Hartree equation (|1.8jl which can be rewritten in terms of the one 
particle Wigner transform as 



d t W t (x,v)+vV x W t (x,v) 



^3 j dydu\({U*p t ){x + £ ^)-{U*p t ){x- E ^))w t {x,u)e™y. 

(3.2) 

Let us denote by W^ n '(t, x, v) the n particle Wigner transform constructed taking tensor products 
of solutions of (|3.2|) . that is 

n 

W( n \t,x,v)=llW t (x J ,v j ). 
3=1 

Moreover we denote by 



H S ^ (t, x, v) = f G% * x G% *„ W%> (t) ) (x, v) and 



H 5 ^ (t, x, v) = ( G£ * a *„ WW ( t ) j ( X , v 



(3.3) 



the Husimi functions associated with the solution W$(t) of ()3.1|) and with W^(t), respectively. 
Here we used, as in Section f2. 11 the notation 



-i \ 3n/2 2 



7T 



<5 2 



(3.4) 



The main result of this paper can be stated as follows. 
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Theorem 3.1. Let U be a radial symmetric real valued potential and assume that there is a constant 
Ki so that 

\\U\\ m = J \Um\Z\ m d£<*T™l (3-5) 
for all m S N. Suppose that, for k < 21ogiV ; 

O( k \w£\Q)-WW(0)) <^sup|0«(x,v)|. (3.6) 



Then, for any fixed I, 6% and 82, we have 



limsup sup 

N— >oo x,v 



uniformly for all t < \{y/l + 1/{7k\) — 1). 
Remarks . 

1) Condition (J3.5|) holds for bounded, real analytic functions U(x). The symmetry condition is 
physically natural. The proof can easily be modified to include non-symmetric potentials as 
well. 

2) It is clear from the proof (see Section 0} that the theorem is still true, with N in (|3.7|) replaced 
by N~ 1+K with an arbitrary small k > 0, if we allow 6%, 62 and £ depend on N as long the 
conditions 

t(N) = o(y/te£N) , [SjiN)}- 2 = o( v / foiiV), j = l,2, (3.8) 

are satisfied. 

3) It is also clear from the proof (see Section that the theorem is still true if we replace the 
accuracy 1/N both in (|3.6|) and (j3.7|l by N~ K with < k < 1. It follows that if the exchange 
term We^ (see Section \l .2(1 of the initial data is of order N~ K , then it remain of the same order 
for all sufficiently small times. 

The proof of this theorem is given in Section 0] below: it is based on a perturbative expansion of 
solutions of the BBGKY Hierarchy. For technical reasons, instead of expanding solutions of ()3.1|) . it 
turns out to be more convenient to work with the Fourier transforms fJi, 17) of the ■ (x, v) (see 
Section [2.31 for the definition of Eq. ()3.1|) is equivalent to the following hierarchy of equations 
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for the marginals nffl'- 



3=1 



-2e 2 Yl [ dqU( q ) S m( £q -^ Vk) ) $ (t, ft, . . . , - g, ■ ■ ■ & + g, ■ ■ ■ , £ n ,rj) (3 . 9) 
- (1 - ne 3 ) I dg?7(g) - sin (^) /^ +1) (t, 6, . . . & - g, . . . £„, <7, *7, 0) 

3.1 Vlasov hiearchy 

The classical Vlasov hierarchy is the semiclassical approximation of the BBGKY hierarchy. It is 
obtained from ()3.1|) by formally setting e — ► and approximate the potential difference by gradient. 
Using 

U (xj + ~ x n +ij ~ U { x 3 ~ ~Y ~ x n+i) = VJ7 (xj - x n+1 )eyj + 0{e 2 ) (3.10) 

and 

we can perform an integration by parts, then integrate out du\dy\ . . . du n dy n to collect delta functions 
S(uj — Vj). Neglecting lower order terms, we obtain formally 

n 

d t W {n) (t; x x , . . . x n , vi, . . . v n ) + Y v o ' v ^ W {n \t- x x ,...,x n ,vi,...,v n ) (3.11) 

n „ 

= vc/ ( x i ~ 2;n+i)V^.^ (n+1) (t;xi, . . . ,x n , x n+1 ,vx, . . . ,v n ,u n+1 )du n+1 dx n+1 

3=1 

for the weak limit 

W (n) (t,x,v) := lim Wt\t,x,v) . (3.12) 

The main result of |14j and ^5] proves that this limit exists, it solves the Vlasov hierarchy (|3. llf) 
and the solution is unique. Therefore = wf where wt satisfies the Vlasov equation (jl.7|) . 

3.2 Other Scalings 

Theorem 13.11 identifies the limit dynamics of the Wigner transform at scale e = iV -1 / 3 . One may 
define the Wigner transform at a different scale v by 

W $U X '> V ) : = JffiS J e~ m - y ^\x + V y -,x-u y -)dy. (3.13) 
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The following lemma shows, however, that under a natural energy conditions, wjp cannot 
converge to a non-trivial function unless v ~ e. Similar statement is true for higher order marginals. 
It justifies our choice of scaling in Section 12.11 in order to derive a dynamics for a nondegenerate 
limiting distribution. 

Lemma 3.2. Let e := iV" 1 / 3 . 

(i) Suppose the kinetic energy of a state \P is comparable with the ground state kinetic energy of 
N fermions in a box of size one, i.e. 

N 

( < C1N 5 / 3 . (3.14) 

Let O £ S(M. 3 xl^) be a Schwarz function with support in {\v\ > A} for some A > 0. Then 

\{0,H NtU )\< [C^Y + O(v\- 2 )\\\O\\ O0 , (3.15) 

with Hn,v := W^ u *x G^-k v G in particular, the one particle Husimi function of ^ at scale 
v vanishes outside of the {v = 0} hyperplane if v -C £. 

(ii) Suppose that the average mean square displacement of ^ is of order one, i.e. 

1 N 
i=i 

Let O G <5(R^ xl^) be a Schwarz function with support in {\v\ < A} for some A > 0. Then 

\{0,H N>V )\ < (const.mOW^-) (3.17) 

with a universal constant. 
Proof, (i) Since \0(x, v)\ < HOIloo^ 2 -^ -2 an d the Husimi function is positive, we obtain 

l^l|oo / 2 tt \ Halloo /' O — . , * , n . IIOI 



\{0,Hn,u)\ < ^ 2 (v ,Hn,v) = ^2 J u G^(v - u)Q v {v)dvdu = J (v + cv)g u (v)dv, 

where q v {v) := J W(x,v)dx is the momentum distribution and c is a universal constant. Since 
/ g u (v)dv = 1 and by (|3~T1|) 

J v 2 g u {v)dv = v 2 Tr(-A) 7 < C x v 2 N 2 ^ = Cxiyjef , 
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we obtain (|3,15|) , 

(ii) We apply the Lieb-Thirring inequality ^2] i n the Fourier space 

r N 

/ g(vf/ 3 dv < (const.) ( xfj ^ (3.18) 

J 3=1 

with a universal constant, where g(v) = N J \ty(v,V2, . . . ,VN)\ 2 dv2 ■ ■ ■ dvjy is the one particle mo- 
mentum density of the antisymmetric function After rescaling we obtain g v (v) = (e / V) 3 g(v / v) , 

c /o ^^^^^^ ^^^^^^ 

hence J gj < (const. )C2(e/^) 2 from (|3.16|) and (|3.18j) . Therefore 

\(0,H NtV )\ < ||0||oo y daxkx(M < X)Hn )V (x,v) 

= ||0||oo y dudu x(|«| < *)Gy*( u - «)a/(«) ^ 

< (constOIIOIUHxd • | < A)|| 5/2 ||G^||i||^|| 5/3 

/ Xe \ 6/5 

< (const.)C 2 ||0||oo(— J 

by Young's inequality. □ 

This lemma shows that the weak limits of Wjp are zero if v e, in particular if the Wigner 
transform is unsealed, v = 1. It may, nevertheless, be reasonable to investigate how well Hartree 
or Hartree-Fock evolutions approximate the true dynamics compared to the actual size of W in a 
different topology. 

Bardos et al. 2 have recently studied the equation (11.10(1 and showed that the Hartree-Fock 
equation approximates the dynamics in the trace norm. In order to study (jl.lOj) one first needs 
to choose the parameter a. Denote by H^, a = —aAjy + (1/N) J2i<j U(x{ — Xj) the Hamiltonian 
corresponding to Ijl.lOj) . and consider an initial state 7jv,o- Here we assume the two body potential 
U to have bounded derivative. Let 7v,t be the time evolution of 7at,o- We are interested in an 
estimate for the mean squared distance between two particles at an arbitrary fixed time t. Define 
the quantities 

u t : = [Tr 7^ [x± - x 2 ) 2 ] 1/2 , 
v t : = [Tr j Nyt (pi - P2) 2 ] 1/2 , 

where pj := —i'V Xi . For typical interacting initial states the mean square distance between the 
particles is of order one, uq = O(l), the kinetic energy per particle is of order TV 2 / 3 , due to Fermi 
statistics, therefore v$ < (const^N 1 / 3 . The next lemma shows that vq is exactly of order iV 1 / 3 for 
any fermionic state localized within an order one distance from the center of mass; in particular there 
cannot be strong velocity correlation between the particles. Then in Lemma 13.41 we show how to use 
the lower bound on vq to give a lower bound on the mean square displacement u 2 . 
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Lemma 3.3. Let 7 be a fermionic N -particle density matrix, Trj = 1, satisfying 

N , N 

: • 

Tr 



7^X>;-*) 2 ]<^ -V: (3-20) 
3=1 3=1 

T/ien 

JV 7 (p! -p 2 ) 2 > (const. )N 2/3 
with a positive constant depending on K . 

Remark. By the symmetry of 7 and a Schwarz inequality 

N , N 



Tr 



l^(x,-X) 2 l=Tr 7 (^E^-^) 



7 AT 



j=i N j=i (3.21) 

< — Tr 7 (xi -X2) , 

so the condition (13.20(1 is satisfied if Tr 7 (x\ — X2) 2 < AT. 

Lemma 3.4. Let C := ||VJ7||oo and let uq, vq be the initial mean squared distance between two 
particles in position and momentum space, respectively. Then for any < t < vq/(8C) 

u 2 > Uq + a 2 v^t 2 — {const. )at(^iQVQ + u$t + avot 2 ^ , (3.22) 

where the constant depends only on C. 

The proofs of these lemmas are deferred to the Appendix. 

According to Lemma 13.31 and the subsequent remark, if the initial inter-particle distance uq is of 
order one, then > {const.)N 1 ^ . In this case Lemma [3 .41 shows that if we want ut to remain of order 
one for t > uniformly as A" —* 00, then we have to assume that a = 0(e) = 0(N^ 1 ^ 3 ). Otherwise 
the interaction between the particles typically vanishes as U(x\ — X2) — > for \x% — X2I — * 00. 

When a = e, we can rewrite the Schrodinger equation Ijl.lOj) as 

i£d t ipN,t = I - e 2 ^2 A :Ej . + ^ U (xj - x k ) j ipN,t- (3.23) 

V j=l j,k J 



This equation is the same as (jl.2|) except the extra e factor in front of the interaction. Since (|1,2[) 
converges to the Vlasov equation, (|3.23|) converges to a free evolution. 

Although some of these conclusions are partly based on initial data considered in Jl], |15| or 
Section 12.21 this behavior is expected for a general reasonable interacting physical system. While 
one may be able to consider some initial data so that the one particle density matrix 7IJI (for the 
dynamics (|1.1U(0 is not given by a free evolution in the A" — ► 00 limit, we do not know if there is a 
natural class of such initial data. 
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4 Proof of the Main Result 



As explained in Section |3] the proof of our main result, Theorem 13.11 is based on a perturbative 
expansion of solutions /xjy (t, £, 77) of the BBGKY hierarchy in the form (|3.9|) . We will compare 
fJ$ (t,£,rj) with tensor products of a solution of the Hartree equation (|3.2|) . which, after Fourier 
transform can be written in the form 

OtMZ, ri)=ti- V^(e, »7) - y dg U(q)^ sin (^) - g, r/)^(g, 0) (4.1) 

with a given initial condition /^o- In the following we will use the notation 

e 

$ ) fori) = ]llH{M (4-2) 

i=l 

for ^"-particle tensor products of a solution /i 4 of (|4.1|l . We remark that global existence, uniqueness 
and regularity of the solution of Q4.1JI have been established in [SJ[7j. 

For any n-particle observable 0( n \£,rj), with £ = (£1, . . . ,£ n ), 77 = (r/i, . . . ,r? n ), we define the 
norms 

„ n 

HO (n) IU= / d^ \0^(t 77)| nd&l + I^D ai 
for q = (ai, . . . , a n ) £ N n . Moreover we use the notation 

(0 (fl) ,M (B) ) = / d^O^fori) M (n) (tr7). 

The following lemma is the main ingredient in the proof of Theorem 13.11 
Lemma 4.1. Assume there exists a constant n\ so that 



(4.3) 



(4.4) 



\U\ 



\U(0m m dC<nTm\ 



for all m £ N. Fix positive integers £, n and suppose that, for all k < (n + I), 



(0,^(0) 



< — llOllo- 



Consider an observable with 

\\0 W \\ a < C^^a 1 \...a e \ , VaeN f , 

then we have 

m TV v m \ 1 — 

where we put Kt = 9K\t(l + 2t)(«i + K2) and assumed that K t < 1. 



^+3 



(4.5) 

(4.6) 
(4.7) 
(4.8) 
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Proof. From (|3.9|) . expanding around the free evolution, we find 

-3 V- t a f a fri ^ 2 ( £ 1 ■ (fa " + (* - s )fe " &)) 



e 3 V [ ds [ dq U(q) - sin 



l<j<k<t 



2 

x fi^'(s,Ci,...^j -q,-..,Ck+q,---, Cn,V + (t- s)t) (4.9) 

^-^Ejf da / d <? 



.0+1) 



2 . feg-Oa + 

sm 



x /% ( s > 6l» • • • , ii - q, ■ ■ ■ , &, q; v + (* - o). 

Next we insert this expansion in the expectation (0^,fj,^S) and we find, moving the free evolution 
from the /i to the observable, 



i<j<k<v 



sm , 

V 2 

X • • • ,£j ~ q, ■ ■ ■ ,Cfe + ?> • • • ! £n! »7) 

- fe3 ) E jf ds / d ^ / d « %)^% »/-(*- s )&- £ sin ( £ ^ 1 



(4.10) 

Now we define the following two operators acting on the observable O^: 

(AO«)(t i) = -e 3 E / dg U{q) 2 - sin ( £g " ( ^~ ?7fc) ) (fc, . . . , & +g, . . . , & - g, • • • , & t,) 

(4.11) 



l<j<k<e 

and 



(BQ M )(f, &+i; %+i) = - E fr(&+i)*fa+i)f sin ' ^ ) 0^(6, • • • , + fr+i, ■ ■ ■ *?)■ 

(4.12) 

Moreover, we denote by (StO™)(£,rj) := 0^(£,rj — tg) the free evolution of the observable O^'. 
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Equation (|4,1U|) can be rewritten as 

(OW /$(*)> = f d£dr, (StOW)(t,T,)f$(P,t,Tl) 

(4.13) 

or, in a more compact form, as 



JU (4.14) 
+ (1 - fe 3 ) / ' ds(BS t - s 0W,^ +1 \s)), 



where we used that the operators A, B and St commute with the complex conjugation (note that, 
since U(x) is symmetric and U(q) is real). Next we iterate this relation n times. We find 

(o«\^(t)) = (s t o(«\^m 

+ dsi ds2 '-- / d Sm (,S SmJ B5 Sm _ 1 _ Sm i?...^_ Sl OW,4 +m) (0)) 

m=l ^0 J JO 

+ d Sl ds 2 ... ds n (BS Sn _ 1 - Sn B...BS t - Sl O^\^N +n \sn)) 
Jo Jo Jo 

n ft fSl pSm—l 

+ V / ds x / ds 2 ... / ds m (AS Sm ^ Sm B...BS^ Sl O^\^ +m - l \s m )) 

m=l JO JO JO 

n f't f'S\ /*S m _i 

-e 3 V(^ + m-l)/ d Sl / d S2 .../ ds m ( J B5 Sm _ 1 _ SmJ B...i?5 t _ Sl OW,4 +m) ( Sm )) . 

m=l - 70 - 70 1/0 

(4.15) 

Using (|2.22|) we have, for any time t and any observable 0^ k \ 

\(0M,$\t))\ < J d$dr, \0Mfor,)\. (4.16) 
So, in order to control the error terms on the last three lines of (|4.15|) we need to estimate the 
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quantities 



| ( H S Sk B S- Sk ftOW) (&, . . . , £n+^; 771, ... , Vn+e ) | d£ dry and 
k=i 

/n—l 
(s Sn AS- Sn [J 5 Sfe B S- Sk S t O®) (6, • • • , Cn+z-i; Vi,'", Vn+i-i) 
k=l 



(4.17) 



d£ dr). 



We begin by i"Q jn . By the definition of the operator B (see ([4.12)0 we have, for general m £ N and 

s£i, 

2 m e 

(S s 5S- s O (ra) )fe,...,( m+ i;!)i,...,i} m+ i) = -^*7(£m+i) sin^|(77j - s£,-)£m+ij 

i=i 

X 5(r/ m+ i - S^m+l) (m) (Cl, • • • , Cj + im+l, ■ ■ ■ , £m5 • • • ,Vj + ■ ■ ■ ,Vm) ■ (4.18) 

Since | sinx| < we obtain the bound 

i+n-l 

£ 

3=1 

H S Sk B S- Sk S t O^)(^,... ,Cn+e-i;Vi,---,Vn+e-i) ■ (4.19) 



/£+n-l 
d£i . . . d£ n+ £dr/i . . . dr/ n+ £_i|f7(^ +n )| \& +n \ ^ |r/j - s„f j 



n-l 



k=l 



Applying equation (|4.18|) once again we find 
Ki >n < / d£i . . .dCn+edm . . . dr] n+ e_i\U (€e+ n )\\&+n\ |J7(& +n _i)| |&+ n _i| <5(^ n+ |_i - s n _i£, t+<? _i) 



>fn-l n+l-2 



n 



31=1 
n-2 



32=1 



| ( II ^ Sfc B S-s k StO^A (£1, ..,£,- 2 + £ n +£-i, ..,^ n+ ^_2; -,r]j 2 + s n _i^ n+ f_i, ..ry n+ ^_ 2 ) 



fe=l 



(4.20) 



After shifting the variables £j 2 — ► £j 2 — £ n -f£-i, % — ^ % ~~ s n-i£n+^-i an d computing the integral 
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(4.21) 



over ?7 n+ £_i (using the delta- function) we get 

Kt, n < Jd^i . ..d£ n+£ dr]i . . .drj n+ t- 2 \U(Ce+n)\\ &+n I |#(&+n-l)| |&+n-l| 

e+n-2 n+l-2 

X ( E l^'i _ Sn ^'il + 2 ( s «-l ~ s n)\€n+i-l\) ( E '/./.! s n ly/V 
ii=i 32=1 

n-2 

x | (J] ^SS^^jfe,..,^^!,...^) . 

fc=i 

After n such iterations we arrive to the estimate 

/n 
d£x . . . d£ n+e drjx... drjt \U(& +k )\ \&+k\ 
k=i 

n £ £+fc-l 

X II (El*» + (*-**)6l+2 E (^-**)fe-|))|0 W «,n)|. (4-22) 

fc=l i=l j=£+l 

Using that |s» — < i for all i,j we get the bound 

/n 
d£i • • • d^ n+ £ dr/i . . . d?# |£/(& +fc )| |£ £+fc | 
fc=l 

n I l+k-l 
X II (E(W+*l6l)+2t E M))\° {i) ^)\- (4.23) 

fc=l i=l j=-£+l 

Let us use the notation 

zi:=E(l»7i|+t|&l); Si == 2*l&+i-il for j = 2, . . . , n . (4.24) 

i=i 

The integrand on the right hand side of equation (|4.22j) is dominated by 

n n k 

( [] • |o w (6, • • • ■ ■ ■ ■ U E>; , ( 425 ) 

fc=l ' k=l j=l 

which, in turn, is bounded by 

n n " i 1 

{X\\u(ik + iWk + i\)^o^(^...^t,m,---,m)\ ■ (Y. n ~ J n x j) . ( 426 ) 

fc=l 3=1 n 
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where we estimated the product by its arithmetic mean in power n. Next we use the binomial 
expansion 

n . , n ( n ~i+ 1 ^.\ aj 

n — j + 1 \™ 



(E^^r- e n 

j=l tti+...+ct»=n j=l J 

and we note that, because of the assumption (|4.7|h we have 

£>I+*|6|J d ^ dr ? 

<Cg(l + t) a «5a! E i<q«(l + t)«*/£j^±^. (4.27) 
This, together with the assumption (|4.5p . implies that 

oi+...+a„=n ' ' j=2 



< eg ^((i + o«i + 2tK 2 r { ^±^ n 



.7 + 1 

3=2 



2 



< c ^^±^ <- V + «2)»(i + a)» (t^),) 2 < »' ( n ^) cfor 1 [9«l(«l + *a)(i + »)]" . 

(4.28) 

Analogously we can bound M^ n . Using the definition of A we find 

M lin < e 3 fdq \U(q)\ \q\ ^ \(rjj - s n &) - (rj k - s n £ k )\ 

j<k 

n-l 

fj S Sr B S- Sr StO®) (&,..., r/i, . . . , »w +n _i) | d£ dr? , (4.29) 



r=l 
and since 

i+n-l 

^2 1^' ~ ~ ~ s n£k)\ <{i + n-2) ^ \Vj ~ Sn£j\ , 

j<k j=l 

we get the bound 

M e , n <e 3 (e + n-2)K e , n . (4.30) 
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Inserting (|4.28|) and (|4,3U|) in (|4,15|) and performing the integration over the s variables, we find 

(ou$®)-{(s t ou$m+ 

+ dsi / ds 2 ... ds m (S Sm BS Sm _ 1 - Sm B...BS^ Sl O^\^ +m \0))\ 

m=1 Jo Jo Jo > 

<{ n+ A < + 2^C V f> + m) ( m + £ ) «f , (4-31) 

\ / m=l 

where we introduced «t := 9tei(/«i + «2)(1 + 2i). Next we want to compare (0^\fi^(t)) with 
(0( £ \(iP), where was defined in (|4.2|1 . Using that fM{£iV) is a solution of the 1-particle Hartree 
Equation (|4.1|) we find that /i^ satisfies the following hierarchy of equation: 



J dqU(q)~ sin (^) /if +1) . . . , & - ?, . . . , £„, g; tj, 0). 



(4.32) 



One can then expand the expectation (O^ , jlf) in a series, exactly as we did for (C>W,/4(t)). Clearly 
one finds 



+ J2 dsi / ds 2 

/O •/ 



1=1 



^m — 1 



d Sm (S Sm BS Sm ^. Sm B . . . BSt- tl O®, +m) > (4.33) 

+ dsj ds 2 ... ds n {BS Sn _ 1 „ Sn B...BS t ^ Sl O^\^). 
Jo Jo Jo 

The error term on the last line can be bounded as before (equation (|2,22|) holds with fj$ replaced 
by JjfJ^ as well). We have 



<0W,^>-{<ft0W,/iW>+ 



m=l 



o Jo 



Si 



n-1 t 

+ Yl / dsi / ds2 • • • 



S-m — 1 



ds m (S Sm BS Sm _^ Sm B . . . BS t ^ Sl O^,^ +m) 

n + i 



< 



! )C V<- (4-34) 
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Combining (jOTj) and ((OH) we find 

<0<'>,(„<$>( t )-A?>)> - Jn + t 
+ |{s 1 o('),(4>(0)-4' , )> 

+ Yl / dsi / ds2 • • • 



\ / m=l ^ ' 



m=l 



(S Sm BS Sm _^ Sm B . . . BS*- n O® (/4 +m) (0) - rf +m) )) 



(4.35) 



Using the assumption Q4.6JI and equation (|4.28|) to bound ||S' Sm i?5 ; j m _ 1 _ Sm S . . . BSt~ Sl O^\\o we find 



Using that 
and that 



< 2 1 ' n + ^ C^r 1 K * + Cfc 3 + 3e 3 C V £ (£ + m) ( m + £ ) k)" . ( 1 .3(5 ) 

m=l 



n 



E«+™)( ra + < )«r s(<+2) '*'£( ra+ m +a ) 

m=l ^ ' m=0 



m + < + a V = <* + *)•«. fT-V)**. (4.37) 



1 — «i. 



the claim of Lemma 14 . 1 1 follows . 



□ 



In order to apply this lemma to prove Theorem I3.1[ we need to estimate the ot-norm of some 
product of Gaussian functions in the and in the 77-space. This is the aim of the following lemma. 

Lemma 4.2. For £ = (£1, . . . , rj = (771, . . . , 77^) we set 



F^ S2 (tv)--= e-^re-^r. (4.38) 
Then there exist universal constants C\ and C2 such that for arbitrary k > 



ipW 11 < ( c i V rW 2 ^ K \<*\ n 1 « 1 



(4.39) 



where 5 1 := S 1 + S 2 , and \ac\ = a± + • • • + ag. 
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Proof. We have 

= / dddr, \Fg] h (tv)\ Um + NP = / d*dn e- 6l f e~*£ {[(^ + \ Vj \ 

j=l i=i 

< ]]2 a >{ j d&e~— |^ y d % e-— + y d^e-— y d^e- — 1^1^} 



(£±Y r(!!i±l 



(4.40) 



for a universal constant Ci. Here we put 5 = (5 1 1 + 5 2 ) 1 - Simple estimate shows that 



thus 



where Ci,Z^2 are universal constants. Elementary calculation shows that 

j-, \ a-j/2 „ „ 



5 2 aj 



< C 2 1/l ° ^ (4.43) 



for a sufficiently large universal constant C%. □ 

We are now ready to proceed with the proof of our main result, Theorem 13. II 

Proof of Theorem I.V.il First we note that the assumption QM.fijl is equivalent to the assumption (|4.fi|) 
in Lemma |4.1l after taking Fourier transform. On the other hand, with the notation 5W^'(t) := 
W$(t)-W®(t) we have 

(fl** - flg*) (t,x, v) = J dx'dv' G^(x - x')G^(v - V)8W^(t,x', v') 



— J I dt&rie^+^e-^- e-^S^Xt,^). (4.44) 



In the following we use the notation 

M ^ ^ ._ f 1 V" i(x-£+v.„) / ^ W 



27 



From Lemma 14.21 we find, for an arbitrary K2 > 

C 



1^*11-^(5^) C'^^W ■■■<**<■ (4.45) 

where the constants C\ and C2 are universal. From Lemma 14. II and from Eq. H4.44|) it follows that 

Ci 



TT&1,82 _ Tj5l,&2 

n l,N n t,N 



(t,x,v) 



< 2k, 



+ 



Ci 

N \2ir5\5l 



' c t/(s 2 4) ( 2Kf y 



a 



1 — Kt 



e+3 



for any K2 > and n < 2 log N — I. Here, as in Lemma 14.11 we use the notation Kt = 9ki (k± + 
K2)t(l + 2t). Since t < \{\Jl + 1/(7k'() — 1), we can fix K2 > such that 2K t < e~ x . Then choosing 
n = log N we find 



£,N n LN 



(t,x, v) 



< 



N 



where Cg^^ is independent of N. Thus, for any fixed £,81,82 we get 



lim sup sup 

N—*oo x v 



TT&1,&2 _ TT&1,&2 

n LN LN 



(t,X, v) 



• X < Ci t Si,5 2 - 



(4.46) 

(4.47) 
□ 



A Proof of Lemma 13.31 and 13.41 

Proof of Lemma VJ.'JX We can restrict ourselves to pure states. Let ^ be a normalized fermionic 
wavefunction. For any l£K 3 define 

• • • , VN-i) ■= ®[yi + X, y 2 + X, . . . , y N ^ + X, X - {y x + . . . + y N - X )) , 

where X := X/N. Clearly is antisymmetric and J*||^x|| 2 dX = 1. By the Lieb-Thirring 
inequality in the Fourier space (|H.18j) 

„ N-l 

/ Qx (v) 5 / 3 dv < (const.)\\^ x \\ 4/3 (^x, ( £ yfj^x) , (A.48) 

where Qx '■= Q^ x 1S ^ ne momentum distribution of the one-particle marginal of *$>x with the nor- 
malization j Qx = (X — l)\\^x\\ 2 (see the proof of Lemma 13.2(1 . Simple calculation shows that 

/ dX (y Xi ( £ y-f) 9 X ) = -^(*, 5>j - A) 2 v&) . (A.49) 
7 i=i ' j=i 



2S 



For an arbitrary p(v) G L 1 (R 3 ) n L 5//3 (M 3 ) we have 



dvv 2 p(v) < jo I dvv 2 p(v), 
'\v\>e J\v\>l t 



f dvp(v) < 4 / 
J\v\>t 1 J\v\. 

/ \3/5 (A.50) 

/ dv P (v) < e 3 U / dV /3 M < ^ 5 ||p|| 5/3 . 

This implies that 

J dvp(v) < iy d V « 2 pW +^ 6/5 |IpII 5 /3- (A.51) 

Optimizing with respect to ■/? we easily obtain J u 2 ^(w)df > (const. )||^||^ 3 /||^||^2 with a positive 
constant. Applying this inequality for qx, using ()A.48|) and the normalization ||f?x||l = (A — 
l)||^x|| 2 , we have 

f 2 , . , (amst.)(N- l) 8 / 3 ||^x|| 4 

Integrating X, using a Schwarz inequality and (|A.49j) we obtain 

f\\*x\\ 2 dX) 



2 



v 2 ^(w)dfdX > (const.)iV 8/3 W ' — * > (const.)N 5/3 



if iV > 2. Finally we conclude by the identity 

JV-l 



/ / v 2 Q X (v)dvdX= ! dx(y x ,(Y,(Pj-PN) 2 )*x) = (N-l)(%(p 1 -p 2 ) 2 V) , 

J J J j=1 

where we again used the symmetry of ^. □ 

Proof of Lemma \S.4\ First we want to prove that vt remains of order iV 1 / 3 for all finite times. 
To this end we compute 

[iH Nj(X , (pi -P2) 2 } = -jj: ^2 ((pi -P2) • (W(xi - x m ) - VU(x 2 - x m )) 

m>3 

+ (VE7 (xi - x m ) - VU(x 2 - x m )) ■ (pi - pa)) 

2 / 

~iV \S Pl ~ P2 ^ ' VC/ ( x i ~ x 2) + VU{x\ - x 2 ) ■ (pi -p 2 ) 
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which implies, using C = ||VZ7||oo, and applying the Schwarz inequality, that 

\8 t v 2 t \ = |Tr {j N , t [iH N>a , (pi -p 2 ) 2 ]) \ < 8Cv t . 
Integrating the last equation we obtain 

v - 4Ct <v t <v + 4Ct (A.52) 
for alH > 0. Next we derive an upper bound for the quantity ut- Here we use 

[iH N , a , (xi - x 2 ) 2 } = 2a({pi - p 2 ) • (xi - x 2 ) + {x\ - x 2 ) ■ (pi - P2)) , 
and from ()A.52|) we find that 

\dtu%\ < 4au t « t < 4 a (v + 4Ct)u t , (A.53) 

hence, for t < vq/8C, 

ut<UQ + 3avQt. (A. 54) 

Finally we want to estimate the quantity ut from below. To this end we compute the second derivative 
of ut using that 

[iH N , a , [iH N ,a, (xi - x 2 ) 2 ]} = 8a 2 (pi -p 2 ) 2 - ^(V?7(xi - x m ) - VU(x 2 - x m )) ■ (xi - x 2 ) 

m>3 

- -j^rVU(xi - X 2 ) ■ (Xl ~ X 2 ). 

Applying the Schwarz inequality, using C = \\S/U\\oo and equations ()A.52|) . I|A.54|) . we find 

d 2 t u\ >8a 2 v 2 -8Cau t 

> 2a 2 Vq — 8C a (uq + 3 a vot) , 

for t < Vq/8C. Integrating this equation twice with the help of (|A.53|) . one easily finds (|3.22|) . □ 
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